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Abstract Recent observations show that the measured 
rates of star formation in the early universe are insuf- 
ficient to produce re-ionization, and therefore, another 
source of ionizing photons is required. In this Letter, 
we examine the possibility that these can be supplied 
by the fast accretion shocks formed around the cores 
of the most massive haloes (10.5 < log M/M & < 12) 
on spatial scales of order 1 kpc. We model the de- 
tailed physics of these fast accretion shocks, and apply 
these to a simple 1-D spherical hydrodynamic accretion 
model for baryonic infall in dark matter halos with an 
Einasto density distribution. The escape of UV pho- 
tons from these halos is delayed by the time taken to 
reach the critical accretion shock velocity for escape of 
UV photons; 220 km s _1 , and by the time it takes for 
these photons to ionize the surrounding baryonic mat- 
ter in the accretion flow. Assuming that in the universe 
at large the baryonic matter tracks the dark matter, we 
can estimate the epoch of re-ionization in the case that 
accretion shocks act alone as the source of UV pho- 
tons. We find that 50% of the volume (and 5-8% of 
the mass) of the universe can be ionized by z ~ 7 — 8. 
The UV production rate has an uncertainty of a fac- 
tor of about 5 due to uncertainties in the cosmological 
parameters controlling the development of large scale 
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structure. Because our mechanism is a steeply rising 
function of decreasing redshift, this uncertainty trans- 
lates to a re-ionization redshift uncertainty of less than 
±0.5. We also find that, even without including the 
UV photon production of stars, re-ionization is essen- 
tially complete by z ~ 5.8. Thus, fast accretion shocks 
can provide an important additional source of ionizing 
photons in the early universe. 

Keywords galaxies: formation - galaxies: halos - cos- 
mology: theory - cosmology: early universe - ultravio- 
let: ISM - shock waves 



1 Introduction 

An understanding of what re-ionizes the Universe at 
the end of the "Dark Ages" has so far eluded us. 
Thus far, two sources of ioniz ing photons have bee n 
suggested; stars and quasars (jLoeb fc Barkanal 120011 1. 
Both of these present d ifficulties. In the case of quasars, 
Mortlock et"all (|2011al ) find very few QSOs at z ~ 6, 
and only one is found above z > 7 (Mortlock, 2011b; 
see also note in proof). More importantly, if there were 
many fainter quasars producing sufficient UV to reion- 
ize the universe at z — 10, then they would overpro- 
duce the present-day unresolved soft X-ray background 
( Dijkstra. Haiman fc Loebll2004h . 

More hopeful is the mechanism of re-ionization by 
massive stars. It is certainly true that star forma- 
tion persists back in time into the early universe. 
The Lyman break techniq ue for detecti ng high red- 



shift galaxies pioneered by ISteidel et alj (|1996l 119991 ) 



has b een successfull y pursu e d out to redshifts o f z ~ 6 



( e. a. iBunker et all d2004h : lYan fc Windhorst! (|2004h 



Bouwens et all (|2006l ): 



Yoshida et al 



(2006)). 



With the 

advent of the Wide Field Camera 3 (WFC3) of the Hub- 
ble Space Telescope this technique has been pushed out 
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to redshifts 2 = 8-9 ([Bunker et al.ll2010l) . In other 
words, we are now directly able to measure star forma- 
tion rates during the epoch of re-ionization. 

Earlier it was thought that the epoch of re-ionization 
occurred at about z rP \ nu ~ 6 — 8. However , the more re- 
cent WMAP results (jrlinshaw et al.ll2009T) place the re- 
ioniz ation as early as z re ion = 10.6 ± 1.4 ([Komatsu et al 
20111) . This causes problems because the measured star 
formation density is found to decrease too quickly with 
increasing redshift. This problem is emphasized by 
Bunker et ail (|2010h . who found that that the Lyman 
continuum photon density arising from star formation 
is in gen eral insufficient t t o reio nize the Universe. More 
recently, iBouwens et al.1 (j 2 1 lh found specifically that 
the ultraviolet flux available from galaxies at z « 10 is 
only ~ 12+io% °f what is needed for stars to be the 
rcionizing source. 

Given the lack of sufficient UV stellar ionizing ra- 
diation - even neglecting the issue of what fraction of 
this escapes from the concentrations of matter in which 
star formation is occurring, we need to consider whether 
there might be another source of ionization which has 
so far been neglected. An obvious candidate is to con- 
vert the mechanical energy of infall associated with the 
growth of structure in the early universe into UV pho- 
tons. An early attempt to est imate UV photon produc- 
tion b y this process was by Sasaki. Takahara fc Sutol 
|l993f ). but their physical model of the UV photon 
production was incomplete, and their mechanism relied 
upon (now ruled-out) isocurvature models, which pro- 
vide much more small-scale structure than the concor- 
da nce ACDM used tod ay. 



Miniati et al.l (|2004l ) considered the production of a 



thermal UV background radiation by turbulence and 
shocks associated with cosmic structure formation, and 
concluded that thermal photons alone are enough to 
produce and sustain He II reionization at z ~ 6. How- 
ever, unlike the scenario considered here, the energy re- 
lease was in the bulk of the baryonic gas rather than be- 
ing concentrated in the highly unusual conditions that 
exist at the cores of the most massive agglomerations 
of mass in the early universe. 

In this paper, we demonstrate that infall of baryonic 
matter into the dense cores of (very rare) massive haloes 
in the early universe creates fast, radiative accretion 
shocks. These generate ionizing radiation in the cool- 
ing zone which can ionize the accreting gas stream and 
escape into the inter-halo space beyond. The theory 
of such "auto-ionizing" fast shocks has previously been 
developed in the context of plasmas of normal chemical 
composition ( Dopita fc Sutherland! 1995 . 1996 ) and ap- 
plied to the study of emission regions a ssociated with 
powerful radio jets and Seyfert galaxies (jBicknell et al 



Il997t lAllen et al.lll998t iBicknell et al.ll2000D . Here, we 
show that fast, "auto-ionizing" accretion shocks are 
very efficient emitters of UV photons. We examine how 
well these shocks, acting alone, could ionize the early 
universe. 

We do not doubt that there may also exist a stellar 
UV continuum to assist the re-ionization, but here we 
will demonstrate that shocks produce a very significant 
number of additional photons at high redshift. In Sec- 
tion 2, we will discuss the radiative properties of these 
fast accretion shocks. In Section 3, we examine the cir- 
cumstances under which such shocks form in the early 
Universe, and in Section 4, we examine the escape of 
these photons from the parent halo. In Section 5 we 
apply this theory to estimating the epoch of accretion 
shock driven re-ionization, and in Section 6 we compare 
this source of photons to the known star formation rate 
at z ~ 8. 



2 Auto-Ionizing Fast Accretion Shocks 



The theory o f fast shocks deve l oped by Dopita fc Sutherland^ 
( 19951 Il996f) and I Allen etaH (|2008l ) cannot be applied 
to the accretion shocks formed in the early universe, 
because the infalling gas is nearly pristine, with a very 
low abundance of heavy elements. Therefore we have 
run a series of fast shock models using the MAPPINGS 
3s code, an updated version of t he co de which was de- 
scribed in[Sutherlan^&^02ita ( 1993h . We have varied 
the shock velocity, v, between 100 and 500 km s _1 . We 
assume the shock to be running into a medium which 
is in collisional ionization equilibrium at 30000 K. The 
abundances adopted are one thousandt h of solar where 



the solar abundances are taken from Grevesse et al 



(2010). These models are assumed to be free of mag- 
netic fields. 

The optical and UV spectrum of two of these mod- 
els is shown in Figure [TJ The spectra are dominated 
by the hydrogen and helium (He I and He II) reso- 
nance and recombination lines, with bound-free and 
two-photon continua superimposed on an underlying 
electron brehmsstrahlung spectrum. These shock spec- 
tra are softer than those computed for plasmas with 
more normal abundance because the channels for cool- 
ing the plasma through the collisonal excitation of 
heavy element resonance lines and continua are absent. 
In particular, these shocks produce many He II pho- 
tons which efficiently ionize H I and He I. The relatively 
'soft' UV spectrum produced in a fast shock with pri- 
mordial composition couples better to the opacity of 
the un-ionized gas, and is therefore more effective in 
producing ionization in the surrounding medium. 
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Fig. 1. — The optical, UV and soft X-ray spectrum 
of fast shocks with velocities of 300 km s~ 1 and 500 
km s _1 . The spectra show the resonance lines of hy- 
drogen and helium, and the 2-photon and bound-free 
continua from these species, superimposed upon a ther- 
mal electron Brehmsstrahlung continuum. Such shocks 
might be detectable in their Lyman-a and confirmed 
by the presence of Hell 1640A emission. 

The cooling timescale of the hot plasma in the shock 
velocity range 200-400 kms -1 is ~ 4 x 10 5 /n yr, where 
n is the hydrogen density in the pre-shock material. 
For the accretion flows we are considering, the density 
at the shock is estimated to be in the range 1 < n < 
10 cm -3 , the cooling timescale is short compared with 
the accretion timescale (~ 5 x 10 7 yr). Furthermore 
the cooling length, ~ 20/n pc, is small compared to 
the physical scale of either the accretion flow, or the 
turbulent disk of shocked gas formed. We can therefore 
treat both the shock and its UV photon emitting zone 
as a thin shell. 

In a plane-parallel shock, 50% of the UV photons go 
to generate a photoionized tail in the cooling zone, while 
the other 50% escape upstream into the accretion flow. 
The number of upstream ionizing photons ^>(u) up pro- 
duced per hydrogen atom processed through the shock 
can be adequately represented from our radiative shock 
models as: 



4>(v)u P = 0.2vf 00 ; 100 < v < 280km/s 
<f)(v) up = 10w| 00 ; 280 < v < 400km/s 
(f>{v)u P = Wvj 0Q ; v > 400km/s 



(1) 



where viqo and V400 are the shock velocities normalized 
to 100 and 400 km s -1 , respectively. 



While 4>(v) up < 1, the upstream photons are ab- 
sorbed in a thin photoionized precursor which is sta- 
tionary with respect to the shock. When <j>(v) up > 1, 
the velocity of the ionization front in the accretion flow 
exceeds the shock velocity and the ionization front is 
free to propagate upstream. Since the density of the ac- 
cretion flow falls with radius from the accretion shock, 
the ionization front will accelerate with time, rapidly 
ionizing the accretion flow and escaping into the low 
density medium. Our accretion shock mechanism for 
re-ionization has the attractive feature that the source 
of the ionization (at the accretion shock) is located 
on the surface of the largest baryonic mass concen- 
trations in the early Universe. The escape of these 
ionizing photons is therefore much more efficient than 
for ionizing photons produced by deeply embedded re- 
gions of star formation. From Equation [T] the critical 
shock velocity for pre-ionization of the accretion flow is 
v w 220km s -1 . This is appreciably larger than the cor- 
responding condition for plasmas of solar composition 
(~ 175 km s _1 ), because the heavy element line and 
continuum cooling of highly ionized species more read- 
ily produce ionizing photons in the hotter zones of the 
shock than does gas of primordial composition. These 
fast shocks produce a radiation field with much higher 
ionization parameter in the precursor gas than would 
stars. The ionization parameters, 0.001 < U < 0.02, 
are comparable with those produced by AGN, although 
their soft X-ray fluxes are much less. 



3 Making UV photons in massive halos 

In our model, the fast accretion shocks needed for re- 
ionization can only be formed around the boundaries of 
the most massive collapsed agglomerations of baryonic 
matter in the early Universe. These mass concentra- 
tions must in turn be located at the cores of the most 
massive halos that can exist at a given epoch. Given 
that there is little star formation, there is also little or 
no feedback to prevent the collection of baryonic matter 
in the core. This (nearly radial) infall of baryonic gas 
is also found in sophis ticated simulation s such as those 
of the OWLS project (|Sales et alJlioich . 

In the standard cold accretion model, the flow re- 
mains cool as long as the compressional heating rate 
and the heating rate by turbulence is matched by the 
cooling rate, locally. In a large fraction of the Universe 
such cold accretion flows have their inflow velocity lim- 
ited by the compressional heating and the need to radi- 
ate this heating away. However, in massive halos, pres- 
sure balance is lost when the cooling timescale becomes 
shorter than the accretion timescale. In 'cold' flows the 
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gas is near 10 4 K and the cooling timescale is of order the 
recombination timescale, r rcc ~ 10 5 /n yr, where n is the 
hydrogen density (cm -3 ). Since the free-fall timescale 
is of order 10 8 yr, the gas loses hydrostatic support and 
free-falls wherever the density exceeds~ 10~ 3 cm" 3 . In 
the halos we are to consider, this occurs at about 6 scale 
radii or 50 kpc. Within this radius, the gas is effec- 
tively pressure-free and in free-fall, achieving velocities 
of several hundred km s _1 by the time it reaches the 
halo core. 

In this picture of rapid accretion onto a dense bary- 
onic core in the dark matter halo, the production of 
ionizing photons, N, is given by: 



M 



(2) 



where M. is the mass accretion rate, p is the mean 
atomic weight of the infalling gas, mjj is the mass of 
a hydrogen atom. For each halo, we therefore need to 
compute both M. and the velocity of accretion onto the 
baryonic core, which by use of equation [T] will then give 
us 4>{v) up . 

For the halos, w e use the Einsato d ensity profile 
( Navarro et~alll2004 ISpringel et "all l2008t ) for the DM 

density, in terms of a scaled radius x = r/r s : 



p(x) 



p-2 exp 



V-i) 



(3) 



In this equation, the logarithmic derivative is taken to 
be equal to —2 at the scale radius, so that p(l) — p_2- 
The parameter a has a small variation about a = 0.18, 
which is the value adopted here, and p_2 is chosen to 
give the halo mean density inside the virial radius, phaio, 
once the physical virial radius and halo over-density are 
known. 

The overall scaling is obtained by requiring phaio to 
be related to the critical density of the Universe at the 
redshift considered, 



OTTLr 



(4) 



Assuming the collapse of a spherical top hat perturba- 
tion, the ratio of the mean halo density to the criti- 
cal density of the Unive rse, A c = ph&io/Pc, is given by 
fervan k Norman|[l99i ); 



A c = 18tt 2 + 82/(2) - 39/(z) 2 , 



where 
/(*) = 



Mi + zf 



- 1 . 



(5) 



(G) 



We take Q M = 0.2715 and Q A = 0.728 from the 
WMAP+BAO+iJo maximu m likelihood value i n the 
WMAP seven year analysis (jKomatsu et al 201 lh . 

The value of the virial radius cc V j r in units of the scal e 
radius is obtained using a fit by iDuffv et alj (|2008f) . 
(termed as concentration in that work) with the desired 
virial mass for the appropriate redshift z, and cosmol- 
ogy (here h = 0.704) : 



x vir = 8.82(1 + z) 



-0.87 



M v 



2.0 x 10 12 h- 



-0.106 



(7) 



Once the virial radius is known in terms of the scale 
radius, the physical scale is chosen so that the mass 
of the halo averaged over the virial volume gives the 
desired mean halo density. The gravitational potential 
is then calculated integrating $(r) = GM(r)/r. An 
analytical form for this potential can be readily derived. 

As an approximation of the complex physics involved 
in the collapse process, we assume that the baryonic 
matter follows the dark matter until virialization of the 
dark matter halo, at which point the baryonic mat- 
ter begins to fall freely towards the centre of the halo. 
Thus, in our initial configuration, the baryonic mass 
starts off from rest and with the same radial density 
profile as the dark matter halo. The baryonic matter is 
therefore already bound within the potential of the dark 
matter halo, and therefore any kinetic energy of infall 
to the initial radius during the assembly of the halo 
is effective assumed to have been already dissipated - 
presumably through the low velocity shocks associated 
with the halo assembly in the cold accretion scenario. 

From this (somewhat uncertain) approximation to 
the initial state, we have followed the subsequent accre- 
tion in 1-D spherical coordinates us ing the Lagrangia n 
hydrodynamic code, Fyris Alpha ( Sutherland 2010l ). 
This code includes all relevant gas physics and gravi- 
tational fields. By using 5 nested domains, a final res- 
olution of 0.3 pc is achieved in this simulation. This 
should be compared with the ~ 0.5 kpc spatial resolu- 
tion achieved in the best smoothed particle hydrody- 
namics (SPH) cosmological codes. 

As expected, within the critical radius at which pres- 
sure support is lost due to fast cooling (within about 
R ~ 50 kpc as estimated above) the baryonic matter 
effectively free-falls towards the centre, since it has no 
thermal pressure support, and very little rotation. In 
the vicinity of the nucleus, the infalling gas is shocked, 
and collects into a small, dense, rotationally-supported 
disk whose size is determined by the torquing of the 
parent halo. The vertical scale height of this disk is de- 
termined by thermal and turbulent support. Its pres- 
sure is given by the ram pressure of the accretion flow, 
since the radiative shocks can be treated as isothermal. 
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Here we assume that the disk has a rotational velocity 
similar to the free-fall velocity in the dark matter halo. 
Thus the depth of the central potential (baryonic plus 
dark matter) is approximately doubled, and the veloc- 
ity associated with the accretion shock is w \^2x the 
free-fall velocity in the dark matter halo. 

For each (non-rotating) spherical halo, the time de- 
pendence of the accreted mass, the free-fall accretion 
rates, and the free-fall velocity at r ~ l.Okpc were cal- 
culated. Using equation [2] we computed the rate of 
production of ionizing photons launched into the accre- 
tion flow. In figure [2] we show a representative halo 
with mass 10 11 Mq within the virial radius. The peak 
mass accretion rate is reached at 50-60 Myr, and cor- 
responds to 205M Q yr~ 1 , when the accretion velocity is 
320 — 350km s _1 . Because we do not properly model 
the initial formation of the halo, it is not clear whether 
photons start to escape the shock contemporaneously 
with this formation, or whether we have to wait an ad- 
ditional ~ 50Myr as indicated in figure O We therefore 
compute both cases. 

The rate of production of ionized photons contin- 
ues to increase with time, as the accretion velocities 
continue to rise. For more massive halos Ai oc M, and 
since v oc M 1 / 2 , equation[5]along with equation[T] shows 
us that N oc M 2 . For less massive halos, the rate of 
production of UV photons falls steeply, due to both the 
steeper dependence of N on velocity, and the fact that 
v(t) < 220 km s _1 for an extended period of time from 
the start of the collapse. 
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Fig. 2. — For a halo with mass lO n M0 within the virial 
radius, we show the mass accretion rate, mechanical 
luminosity through the accretion shock, and the rate 
of production of ionizing photons as a function of time 
after the halo first forms. Note (at late times) how a 
decreasing mass accretion rate, M, is compensated by 
an increase in the accretion velocity so that the rate 
of production of ionizing photons remains essentially 
invariant with time. However, the number of ionizing 
photons produced per baryon accreted is only of order 
0.1% of what would be produced if these baryons were 
all formed into stars. 



4 The escape of photons from the halos 

It is not enough for the accretion shocks to produce 
an excess of UV photons to provide re-ionization. The 
photons have to first ionize the gas in the halo, and 
even when this has been done, we still have to account 
for the number of photons being used per unit time in 
maintaining the ionization of the halo. In our numer- 
ical models, we find that, over a wide range in radius, 
the number density n(r) of infalling gas in our Einsato 
density profile dark matter halos can be approximated 
by a power law n(r) oc r^ 1 - 83 . Since the number of re- 
combinations per unit volume is oc n(r) 2 , the number 
of photons needed to maintain ionization in a spher- 
ical shell of unit thickness in the accretion flow falls 
off oc r~ 166 . Thus, the photons required to maintain 
ionization rapidly become unimportant, and once the 
ionization front moves away from the vicinity of the 
accretion shock, it rapidly accelerates outwards. 

The ionizing photons must still ionize all of the bary- 
onic matter in the halo before they can escape into 



the universe at large. The number of photons required 
is considerable. For a halo with mass 10 11 M Q within 
the virial radius, the total baryonic mass to infinity is 
1.4 x 10 11 M Q . We compute that, once the UV photons 
start escaping from the vicinity of the shocks (~ 50 Myr 
after the formation of the halo), it takes a further ~ 60 
Myr to complete the ionization of the baryonic halo. 
For more massive halos, N oc M 2 so that the timescale 
taken to ionize the halo falls oc M~ x . The steep mass 
dependence of the photon production coupled with the 
decreasing ionization latency time with mass allows the 
more massive halos to successfully compete for the net 
UV photon production, despite their much lower space 
densities. 

Detailed hydrodynamic modeling o f the baryonic gas 
( Greif et alj|2008t iGoerdt et allboioh shows that it is 
accreted in streams, rather than isotropically. This is 
the so-called "cold accretion" model. In this, the early 
escape of photons is favored. This is because the spe- 
cific photon flux scales as the mass accretion rate per 
unit area, so more photons per unit area are produced 
in denser accretion streams. These photons can then 
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easily escape between these accretion streams, result- 
ing in more prompt re-ionization. 




10.0 9.0 8.0 7.0 

Red shift (z) 

Fig. 3. — The fraction of the volume, and of the to- 
tal mass, of the universe that is ionized by accretion 
shock-produced UV photons as a function of redshift,. 
We use the new high-re solution cosmological calcula- 
tion o f Kobayashi, (see Kobavashi. Springel fc White] 
( 20071) 'or the description of an earlier version of this 
code). The ionization model used is a photon diffusion 
approximation in which the ionization fronts are as- 
sumed to move rapidly to ionize the least-dense phases 
first. The number of ionizing photons needed to main- 
tain the ionization in the ionized volume are explicitly 
included in this calculation. We have computed two 
cases - that in which the shocks photon production is 
contemporaneous with the formation of the halo (left- 
hand curves) and the case in which this photon produc- 
tion is delayed by 50 Myr (right-hand curve). Accord- 
ing to our model, 50% of the volume (5-8% by mass) 
of the universe is ionized by z ~ 8, and ionization is 
essentially complete by z = 5.8. 



5 Reionizing the Universe 

The rate of production of ionizing photons be co-moving 
Mpc 3 can be calculated using the density of halos as a 
function of mass and redshift given by, for example, 
the sim ulations of Mo fc White! ( 20021 ) or Jenkins et al. 
( 200lh . EquationElgives the photon production rate for 
a halo of a given mass. We have computed two cases 
- that in which the UV photon production in shocks is 
contemporaneous with the formation of the halo and 



the case in which this photon production is delayed by 
50 Myr, related to the free-fall collapse time. 

The photon production rate is summed o v er the 
space density of halos given by iMo fc White! d2002h 
to give the production rate of photons per co-moving 
Mpc 3 . We then compute the time dependent ioniza- 
tion of the universe using a photon diffusion approx- 
imation. In this, the ionization fronts move rapidly 
to ionize the least-dense phases first, and the denser 
HI in the halos is ionized more slowly and at later 
times. This provides a reasonable ap proximation to 



this complex process (|Iliev et al.l 120091 ) . The photons 



needed to maintain the ionization in the ionized vol- 
ume are included in the calculation. The appropriate 
densities have been computed for us using a more re- 
cent (highe r resolution) version of the cosmologic al sim- 
ulation by iKobavashi. Springel fc White! (j2007h . The 
result is shown in Figure O We find that with the 
Mo fc White! ([2002) cosmological parameters, ~ 50% 



by volume (and 5-8% by mass) of the universe is ion- 
ized by z < 8, and that re-ionization is complete by 
z > 5.8. Note that the effect of more prompt escape of 
UV photons would be to markedly increase the degree 
of ionization at early times. At a redshift of z = 10 
the photon density is increased by more than a factor 
of three by this prompt escape. 

Note that the space density of mas sive halos given by 
the equations oflJenkins et al.l (|200ll ) and the models of 
Mo fc White! |2002) assume the cosmological parameter 



<7g. This is related to the mean mass density fluctua- 
tions on scales of 8/i _1 Mpc today, and is ss 0.9. This 
enters exponentially into the overall collapse fraction. 
There is some uncertainty in this valu e and the more re- 
cent estimates iKomatsu et all (|201 lh suggest cr 8 w 0.8. 
If one assumes this latter value, the fraction of baryons 
in massive halos at high redshift will be reduced by ap- 
proxi mately a facto r of 4.5 compared to those estimated 
using Mo fc White ( 2002f ). As a consequence, the rate 
of production of ionizing photons will be reduced by 
a similar factor. However, the effect on the epoch of 
re-ionization is much smaller, as we now show. 

As described above, only the most massive halos are 
important as sources of ionizing photons. We estimate 
that dominant halos for production of UV photon s 
have masses in the range 10 11 " 11 - 5 M Q . The photon 
production by accretion shocks increases very rapidly 
with decreasing redshift, thanks to the rapid increase in 
space density of halos, and the increasing importance 
of the more massive halos. At z ~ 8, it is increas- 
ing (with decreasing redshift) oc (1 + z)~ 10 . At the 
same time the number of photons needed to ionize the 
universe decreases oc (1 + z) 3 . These factors serve to 
"lock in" the epoch of re-ionization into a rather nar- 
row range in z. Furthermore, as a result of this strong 



Re-ionizing the Universe without Stars 



7 



redshi ft dependence, changing from the Mo fc White! 
( 2002 f) cosmological p arameters to the ones favored by 
Komatsu et all ((20111 ) would reduce the onset of ioniza- 
tion in redshift by only 5z « 0.5. 



6 Comparison with star formation rates 

The accretion shocks provide 7 x 10 52 ionizing pho- 



tons s for a halo mass of log M/M i 



© 



1 1 , and 



as many as 8 x 10 ionizing photons s for a halo 
mass of log M/Mq — 12.0. These fluxes are equiv- 
alent to the ionizing photons produced by 1400 and 
160,000 05 stars, respectively. These numbers may 
seem high, but correspond to the relatively modest 
total star formation rates of 0.07 Mq yr" 1 and 8 
Mq yr _1 , respectively. Allowing for the space den- 
sities of the halos of the appropriate mass, at a red- 
shift of z — 8 the UV photon density provided by 
a halo mass of log M/Mq = 11 is equivalent to 
a star formation density of 1O -3 ' 75 M yr" 1 Mpc -3 . 
For a halo mass of log M/Mq = 11.5, this rises to 
10- 3 - 2 Af Q y^ 1 Mpc~ 3 and by log M/Mq = 12, has 
fallen to 10- 41 M Q yr" 1 Mpc -3 . The total effective 
equivalent star formation density of shocks at z = 7 — 8 
is (1 - 2) x 10~ 3 M Q yr" 1 Mpc" 3 , between half and all 
of the currently measured star formation density at this 



ol trie currently measured, star tormation (tensity at tins 
epoch ((Bunker et alJlioioHBouwens et al.ll201lh . Thus, 



even if accretion shocks do not dominate the UV photon 
budget in the early universe, they remain an important 
contributor to the total. 



7 Conclusions 

We have shown that re-ionization can, in principle, 
be driven by the UV photons produced by fast, ra- 
diative and "auto-ionizing" accretion shocks associated 
with the most massive haloes formed at any epoch in 
the early Universe. These shocks are located at the 
outer edges of the most massive agglomerations of bary- 
onic matter, and therefore the UV radiation produced 
in such shocks can readily escape into inter-galactic 
space. According to this model, the minimum mass 
halo that can produce an appreciable UV flux has a 
mass £ 3 x 10 10 A/q. Such halos are not formed at ap- 
preciable space density until z ~ 12, so the epoch of 
re-ionization cannot start before this. The production 
of ionizing photons increases very rapidly thereafter, 
while the number of photons required to complete re- 
ionization falls. Thus, this process provides a epoch of 
re-ionization which is rather sharply defined. We find 
that 50% of the volume (and 5-8% by mass) of the uni- 
verse is ionized by z < 7 — 8, and that re-ionization is 



complete by z > 5.8. The UV photon density produced 
by these accretion shocks at z ~ 7 — 8 is approximately 
half that produced by the observed star formation at 
this redshift, so even if accretion shocks do not domi- 
nate, they remain an important contributor to the total 
UV photon production budget. In addition, the accre- 
tion shocks are formed around the periphery of the main 
mass agglomeration, and so the escape of the UV pho- 
tons is favored, especially when the gas accretion is in 
cold streams. 

The stars, on the other hand, will be predominantly 
formed within these dense regions where the pressure 
is of order P/k ~ 3 x 10 s cm 3 K. Under these condi- 
tions, individual the H II regi ons around O or B stars 
are only ~ 0.5 pc in diameter (|Dopita et al J 120061 ). In- 
dividually these stars will be ineffective in ionizing the 
universe. Collectively, they could blow out of the dense 
disk gas, and UV photons could then escape. However, 
the timescale for this to occur is appreciably longer than 
the baryon free-fall time and is unlikely to impact upon 
our estimate for the epoch of re-ionization. A doubling 
of the UV radiation density to include a possible stellar 
contribution would push our estimate of the epoch of 
re-ionization redshift back by z ~ 0.6 in redshift. 

To resolve the accretion shocks will take much higher 
resolution than has hitherto been achieved in simula- 
tions. The thickness of the hot shocked gas is only 
~ lOpc. Current cosmological models have typical par- 
ticle sizes of 10 6 M© and a maximum spatial resolution 
~ 500pc. Since the infalling gas has a cooling timescale 
much shorter than the infall timescale, it is not hotter 
than 10 4 K. The post-shock turbulent gas has a similar 
temperature, and so current simulations will make it 
appear that the gas is accreted entirely in "cold mode". 

More precise estimates of the epoch of re-ionization 
will have to await numerical hydrodynamic simulations 
which correctly compute the temporal and mass de- 
pendence of the baryonic fast accretion shocks, which 
explicitly compute the propagation of the ionization 
fronts, which compute the UV photon production rate 
from individual halos, and which correctly model the 
growth of the nuclear disk including disk rotation and 
turbulence. Such models are in preparation, and these 
models will serve to determine to what extent these 
proto-galactic shocks are observable, either through 
their Lyman-a or their Hell 1640A emission lines or 
through the soft thermal X-ray continuum which they 
produce. 
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Note Added in proof: 

The discovery of a massive quasar at a redshift of 
z = 7.085 by Mortlock et al. (2011b) has a number 
of important implications to both this paper, and to 
existing comological models. First, as stated by the 
authors, it establishes that there is an appreciable un- 
ionized fraction at this redshift, in contrast with the 
WMAP result (Komatsu et al. 2011). Second, its spec- 
trum establishes that the broad line region has chemical 
abundances similar to solar, even though the universe 
is only 0.77 Gyr old at this epoch, implying a large 
conversion of baryons into heavy elements in its vicin- 
ity. Third, the estimated black hole mass, 2 x 10 9 M Q 
implies very prompt collapse. If built up at the Edding- 
ton rate from an initial mass of order IOOMq, we would 
need about 3 x 10 s yr to do this! This implies that the 
host halo, M ~ 10 n M o was not only collapsed but 
had formed a super-dense core by z ~ 10, supporting 
the possibility of monolithic collapse of massive halos in 
the early universe. Monolithic collapse, the production 
of massive black holes, as well as possibly concentrated 
regions of star formation suggests that all three sources 
of UV photons may be important for re-ionization, and 
that the current paradigm of massive halo formation is 
likely incomplete 
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